Models of quantum computation and quantum programming languages 



Jarosiaw Adam MISZCZAK 

Institute of Theoretical and Applied Informatics, 
Polish Academy of Sciences, Baltycka 5, 44-100 Gliwice, Poland 

The goal of the presented paper is to provide an introduction to the basic computational mod- 
els used in quantum information theory. We review various models of quantum Turing machine, 
quantum circuits and quantum random access machine (QRAM) along with their classical counter- 
parts. We also provide an introduction to quantum programming languages, which are developed 
using the QRAM model. We review the syntax of several existing quantum programming languages 
and discuss their features and limitations. 
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I. INTRODUCTION 

Computational process must be studied using the fixed 
model of computational device. This paper introduces 
the basic models of computation used in quantum in- 
formation theory. We show how these models are defined 
by extending classical models. 

We start by introducing some basic facts about clas- 
sical and quantum Turing machines. These models help 
to understand how useful quantum computing can be. 
It can be also used to discuss the difference between 
quantum and classical computation. For the sake of com- 
pleteness we also give a brief introduction to the main res- 
ults of quantum complexity theory. Next we introduce 
Boolean circuits and describe the most widely used model 
of quantum computation, namely quantum circuits. We 
focus on this model since many presented facts about 
quantum circuits are used in the following sections. Fi- 
nally we introduce another model which is more suited for 
defining programming languages operating on quantum 
memory — quantum random access machine (QRAM). 

We also describe selected examples of the existing 
quantum programming languages. We start by formu- 
lating the requirements which must be fulfilled by any 
universal quantum programming language. Next we de- 
scribe languages based on imperative paradigm - QCL 
(Quantum Computation Language) and LanQ. We also 
describe recent research efforts focused on implementing 
languages based on functional paradigm and discuss the 
advantages of a language based on this paradigm. As the 
example of functional quantum programming language 
we present cQPL. 

We introduce the syntax and discuss the features of the 
presented languages. We also point out their weaknesses. 
For the sake of completeness a few examples of quantum 
algorithms and protocols are presented. We use these 
examples to introduce the main features of the presented 
languages. 

Note that we will not discuss problems related to the 
physical realisation of the described models. We also do 
not cover the area of quantum error correcting codes, 
which aims to provide methods for dealing with decoher- 
ence in quantum systems. For an introduction to these 
problems and recent progress in this area see e.g. [TJ[2]. 

One should be also aware that the presented overview 



of existing models of quantum computation is biased 
towards the models interesting for the development of 
quantum programming languages. Thus we neglect some 
models which are not directly related to this area (e.g. 
quantum automata or topological quantum computa- 
tion) . 



A. Quantum information theory 

Quantum information theory is a new, fascinating field 
of research which aims to use quantum mechanical de- 
scription of the system to perform computational tasks. 
It is based on quantum physics and classical computer 
science, and its goal is to use the laws of quantum mech- 
anics to develop more powerful algorithms and protocols. 

According to the Moore's Law H] the number of 
transistors on a given chip is doubled every two years 
(see FigureJT]) . Since classical computation has its natural 
limitations in the terms of the size of computing devices, 
it is natural to investigate the behaviour of objects in 
micro scale. 
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Figure f . Illustration of Moore's hypothesis. The number of 
transistors which can be put on a single chip grows exponen- 
tially. The squares represent microprocessors introduced by 
Intel Corporation 4 . The dotted line illustrates the rate of 
growth, with the number of transistors doubling every two 
years. 

Quantum effects cannot be neglected in microscale and 
thus they must be taken into account when designing fu- 
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ture computers. Quantum computation aims not only at 
taking them into account, but also at developing methods 
for controlling them. Quantum algorithms and protocols 
are recipes how one should control quantum system to 
achieve higher efficiency. 

Information processing on quantum computer was first 
mentioned in 1982 by Feynman [S]. This seminal work 
was motivated by the fact that simulation of a quantum 
system on the classical machine requires exponential re- 
sources. Thus, if we could control a physical system at 
the quantum level we should be able to simulate other 
quantum systems using such machines. 

The first quantum protocol was proposed two years 
later by Bennett and Brassard [6]. It gave the first ex- 
ample of the new effects which can be obtained by using 
the rules of quantum theory for information processing. 
In 1991 Ekert described the protocol [7] showing the us- 
age of quantum entanglement [8] in communication the- 
ory. 

Today we know that thanks to the quantum nature 
of photons it is possible to create unconditionally se- 
cure communication links [3] or send information with 
the efficiency unachievable using classical carriers. Dur- 
ing the last few years quantum cryptographic protocols 
have been implemented in real-world systems. Quantum 
key distribution is the most promising application of 
quantum information theory, if one takes practical ap- 
plications PHim m t° account. 

On the other hand we know that the quantum mech- 
anical laws of nature allow us to improve the solution 
of some problems [T2"HT1] , construct games [T5J [IB] and 
random walks [171 118] with new properties. 

Nevertheless, the most spectacular achievements in 
quantum information theory up to the present moment 
are: the quantum algorithm for factoring numbers and 
calculating discrete logarithms over finite field proposed 
in the late nineties by Shor [13] . The quantum al- 
gorithm solves the factorisation problem in polynomial 
time, while the best known probabilistic classical al- 
gorithm runs in time exponential with respect to the size 
of input number. Shor's factorisation algorithm is one of 
the strongest arguments for the conjecture that quantum 
computers can be used to solve in polynomial time prob- 
lems which cannot be solved classically in reasonable (i.e. 
polynomial) time. 

Taking into account research efforts focused on dis- 
covering new quantum algorithms it is surprising that 
for the last ten years no similar results have been ob- 
tained [ini [20] ■ One should note that there is no proof 
that quantum computers can actually solve NP-complete 
problems in polynomial time |21L I22j . This proof could 
be given by quantum algorithms solving in polynomial 
time problems known to be NP-complete such as k- 
colorability. The complexity of quantum computation 
remains poorly understood. We do not have much evid- 
ence how useful quantum computers can be. Still much 
remains to be discovered in the area of the relations 
between quantum complexity classes such as BQP and 



classical complexity classes like NP. 

B. Progress in quantum algorithms 

Due to the slow progress in discovering new quantum 
algorithms novel methods for studying the impact of 
quantum mechanics on algorithmic problems were pro- 
posed. 

The first of these methods aims at applying the rules 
of quantum mechanics to game theory [IB, 23 . Classical 
games are used to model the situation of conflict between 
competing agents. The simplest application of quantum 
games is presented in the form of quantum prisoners di- 
lemma |15j . In this case one can analyse the impact of 
quantum information processing on classical scenarios. 
On the other hand quantum games can be also used to 
analyse typical quantum situations like state estimation 
and cloning [2~i] . 

Quantum walks provide the second promising method 
for developing new quantum algorithms. Quantum 
walks are the counterparts of classical random walks 
[171 [TH] , For example, in [25] the quantum algorithm 
for element distinctness using this method was proposed. 
It requires 0(n 2 / 3 ) queries to determine if the input 
{xi, . . . , x n } consisting of n elements contains two equal 
numbers. Classical algorithm solving this problem re- 
quires 0(n\ogn) queries. The generalisation of this al- 
gorithm, with applications to the problem of subset find- 
ing, was described in [26] . Other application of quantum 
walks include searching algorithms [27] and subset find- 
ing problem. It was also shown that quantum walks can 
be used to perform a universal quantum computation 
[2"51 |2"5] . In [30] the survey of quantum algorithms based 
on quantum walks is presented. More information con- 
cerning recent developments in quantum walks and their 
applications can be found in |31j . 

One should note that the development of quantum al- 
gorithms is still a very lively area of research [20l 132] . 
General introduction to quantum algorithms can be 
found in [33J . The in-depth review of the recent results 
in the area of quantum algorithms for algebraic problems 
can be found in |34| . 

II. COMPUTABILITY 

Classically computation can be described using vari- 
ous models. The choice of the model used depends on 
the particular purpose or problem. Among the most im- 
portant models of computation we can point: 

• Turing Machine introduced in 1936 by Turing 
and used as the main model in complexity the- 
ory [35] . 

• Random Access Machine 36, 37J which is the 
example of register machines; this model cap- 
tures the main features of modern computers and 
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provides a theoretical model for programming lan- 
guages. 

• Boolean circuits [38 defined in terms of logical 
gates and used to compute Boolean functions / : 
{0, l} m h- > {0, 1}™; they are used in complexity the- 
ory to study circuit complexity. 

• Lambda calculus defined by Church 3!) and used 
as the basis for many functional programming lan- 
guages [10]. 

• Universal programming languages which are 
probably the most widely used model of computa- 
tion ju. 

It can be shown that all these models are equivalent 
[551 [35]. I n other words the function which is computable 
using one of these models can be computed using any 
other model. It is quite surprising since Turing machine 
is a very simple model, especially when compared with 
RAM or programming languages. 

In particular the model of a multitape Turing machine 
is regarded as a canonical one. This fact is captured by 
the Church- Turing hypothesis. 

Hypothesis 1 (Church- Turing) Every function 
which would be naturally regarded as computable can be 
computed by a universal Turing machine. 

Although stated as a hypothesis, this thesis is one of 
the fundamental axioms of modern computer science. A 
Universal Turing machine is a machine which is able to 
simulate any other machine. The simplest method for 
constructing such device is to use the model of a Turing 
machine with two tapes (35] . 

Research in quantum information processing is motiv- 
ated by the extended version of Church- Turing thesis for- 
mulated by Deutsch |42) . 

Hypothesis 2 (Church- Turing-Deutsch) Every 
physical process can be simulated by a universal comput- 
ing device. 

In other words this thesis states that if the laws of phys- 
ics are used to construct a Turing machine, this model 
might provide greater computational power when com- 
pared with the classical model. Since the basic laws of 
physics are formulated as quantum mechanics, this im- 
proved version of a Turing machine should be governed 
by the laws of quantum physics. 

In this section we review some of these computational 
models focusing on their quantum counterparts. The dis- 
cussion of quantum programming languages, which are 
based on the quantum random access machines (QRAM), 
is presented in Section |Hl[ 

We start be recalling the basic facts concerning a Tur- 
ing machine. This model allows to establish clear notion 
of computational resources like time and space used dur- 
ing computation. It is also used to define other models 
introduced in this section precisely. 



On the other hand for practical purposes the notion 
of Turing machine is clumsy. Even for simple algorithms 
it requires quite complex description of transition rules. 
Also, programming languages defined using a Turing ma- 
chine [43] . have rather limited set of instructions. Thus 
we use more sophisticated methods like Boolean circuits 
and programming languages based on QRAM model. 



A. Turing machine 

The model of a Turing machine is widely used in clas- 
sical and quantum complexity theory. Despite its simpli- 
city it captures the notion of computability. 

In what follows by alphabet A = {aj., . . . , a n } we mean 
any finite set of characters or digits. Elements of A are 
called letters. Set A k contains all strings of length k 
composed from elements of A. Elements of A k are called 
words and the length of the word w is denoted by \w\. 
The set of all words over A is denoted by A*. Symbol 
e is used to denote an empty word. The complement of 
language L C A* is denoted by L and it is the language 
defined as L = A* — L. 



1. Classical Turing machine 

A Turing machine can operate only using one data 
structure - the string of symbols. Despite its simplicity, 
this model can simulate any algorithm with inconsequen- 
tial loss of efficiency [35]. A Classical Turing machine 
consists of 

• an infinitely long tape containing symbols from the 
finite alphabet A, 

• a head, which is able to read symbols from the tape 
and write them on the tape, 

• memory for storing programme for the machine. 

The programme for a Turing machine is given in terms 
of transition function S. The schematic illustration of a 
Turing machine is presented in Figure [2] 

Formally, the classical deterministic Turing machine is 
defined as follows. 

Definition 1 (Deterministic Turing machine) A 

deterministic Turing machine M over an alphabet A is 
a sixtuple (Q, A,S,qa,q a ,q r ), where 

• Q is the set of internal control states, 

• qo , q a , q r € Q are initial, accepting and rejecting 
states, 

• 5: QxAt-^-QxAx { — 1,0, 1} is a transition 
function i.e. the programme of a machine. 
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By a configuration of machine M we understand a 
triple (qi,x,y), qi G Q, x,y € A*. This describes a 
situation where the machine is in the state qi, the tape 
contains the word xy and the machine starts to scan the 
word y. If x — x' and y — b\y' we can illustrate this 
situation as in Figure [2] 
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Figure 2. Computational step of the Turing machine. Con- 
figuration (qi,x'a,biy') is presented in (a). If the transition 
function is defined such that 5(qi, bi) = (g2, 62, —1) this com- 
putational step leads to configuration (qj, x' , fflfoy') (see (b)). 

The transition from the configuration c\ to the config- 
uration C2 is called a computational step. We write chc' 
if 5 defines the transition from c to d. In this case d is 
called the successor of c. 

A Turing machine can be used to compute values of 
functions or to decide about input words. The com- 
putation of a machine with input w G A* is defined 
as a sequence of configurations Cq , C\ , c% , . . ., such that 
Co = (qi,e,w) and Cj h Cj+i. We say that computation 
halts if some Ci has no successor or for configuration Ci, 
the state of the machine is q a (machine accepts input) or 
q r (machine rejects input). 

The computational power of the Turing machine has its 
limits. Let us define two important classes of languages. 

Definition 2 A set of words L G A* is a recursively enu- 
merable language if there exists a Turing machine accept- 
ing input w iff w G L. 

Definition 3 A set of words L G A* is a recursive lan- 
guage if there exists a Turing machine M such that 

• M accepts w iff w G L, 

• M halts for any input. 

The computational power of the Turing machine is lim- 
ited by the following theorem. 

Theorem 1 There exists a language H which is recurs- 
ively enumerable but not recursive. 



Language H used in the above theorem is defined in 
halting problem [35]. It consists of all words composed 
of words encoding Turing machines and input words for 
these machines, such that a particular machine halts on a 
given word. A universal Turing machine can simulate any 
machine, thus for a given input word encoding machine 
and input for this machine we can easily perform the 
required computation. 

A deterministic Turing machine is used to measure 
time complexity of algorithms. Note that if for some lan- 
guage there exists a Turing machine accepting it, we can 
use this machine as an algorithm for solving this problem. 
Thus we can measure the running time of the algorithm 
by counting the number of computational steps required 
for Turing machine to output the result. 

The time complexity of algorithms can be described 
using the following definition. 

Definition 4 Complexity class TIME(/(n)) consists of 
all languages L such that there exists a deterministic Tur- 
ing machine running in time fin) accepting input w iff 
w G L. 

In particular complexity class P defined as 

P = (jTIME(n fc ), (1) 

k 

captures the intuitive class of problems which can be 
solved easily on a Turing machine. 

2. Nondeterministic and probabilistic computation 

Since one of the main features of quantum computers 
is their ability to operate on the superposition of states 
we can easily extend the classical model of a probab- 
ilistic Turing machine and use it to describe quantum 
computation. Since in general many results in the area 
of algorithms complexity are stated in the terms of a 
nondeterministic Turing machine we start by introducing 
this model. 

Definition 5 (Nondeterministic Turing machine) 

A nondeterministic Turing machine M over an alphabet 
A is a sixtuple {Q,A,8,qQ,q a ,q r ), where 

• Q is the set of internal control states, 

• qo,q a ,qr G Q are initial, accepting and rejecting 
states, 

• ScQxAxQxAx {—1, 0, 1} is a relation. 

The last condition in the definition of a nondetermin- 
istic machine is the reason for its power. It also requires 
to change the definition of acceptance by the machine. 

We say that a nondeterministic Turing machine ac- 
cepts input w if, for some initial configuration (qi,e,w), 
computation leads to configuration (q a , 01,02) for some 
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Figure 3. Schematic illustration of the computational paths 
of a nondeterministic Turing machine [35] ■ Each circle rep- 
resents the configuration of the machine. The machine can 
be in many configurations simultaneously. 

words a\ and a 2 . Thus a nondeterministic machine ac- 
cepts the input if there exists some computational path 
defined by transition relation 8 leading to an accepting 
state q a . 

The model of a nondeterministic Turing machine is 
used to define complexity classes NTIME. 

Definition 6 Complexity class NTIME(/(n)) consists 
of all languages L such that there exists a nondetermin- 
istic Turing machine running in time f(n) accepting in- 
put w iff w G L. 

The most prominent example of these complexity classes 
is NP, which is the union of all NTIME(n fe ), i.e. 

NP = (J NTIME (n k ). (2) 

k 

A nondeterministic Turing machine is used as a theor- 
etical model in complexity theory. However, it is hard to 
imagine how such device operates. One can illustrate the 
computational path of a nondeterministic machine as in 
Figure [3] [3S]. 

Since our aim is to provide the model of a physical 
device we restrict ourselves to more realistic model. We 
can do that by assigning to each element of relation a 
number representing probability. In this case we obtain 
the model of a probabilistic Turing machine. 

Definition 7 (Probabilistic Turing machine) A 

probabilistic Turing machine M over an alphabet A is a 
sixtuple (Q, A,5,qa,q ai q r ), where 

• Q is the set of internal control states, 

• qa , q a , q r G Q are initial, accepting and rejecting 
states, 

• 5 : QxAxQxAx { — 1. 0, 1} n- [0, 1] is a transition 
probability function i.e. 



For a moment we can assume that the probabilities 
of transition used by a probabilistic Turing machine can 
be represented only by rational numbers. We do this to 
avoid problems with machines operating on arbitrary real 
numbers. We will address this problem when extending 
the above definition to the quantum case. 

The time complexity of computation can be measured 
in terms of the number of computational steps of the Tur- 
ing machine required to execute a programme. Among 
important complexity classes we have chosen to point 
out: 

• P - the class of languages for which there exists 
a deterministic Turing machine running in polyno- 
mial time, 

• NP - the class of languages for which there ex- 
ists a nondeterministic Turing machine running in 
polynomial time, 

• RP - the class of languages L for which there exists 
a probabilistic Turing machine M such that: M 
accepts input w with probability at least \ if w G L 
and always rejects w if w £ L, 

• coRP - the class of languages L for which L is in 
RP, 

• ZPP RP n coRP. 

More examples of interesting complexity classes and 
computational problems related to them can be found 
in 01]. 

3. Quantum Turing machine 

A quantum Turing machine was introduced by Deutsch 
[42] . This model is equivalent to a quantum circuit model 
[451 146) . However, it is very inconvenient for describing 
quantum algorithms since the state of a head and the 
state of a tape are described by state vectors. 

A quantum Turing machine consists of 

• Processor: M 2-state observables {rii\i Gz M }- 

• Memory: infinite sequence of 2-state observables 
{TOi|l G Z}. 

• Observable x, which represents the address of the 
current head position. 

The state of the machine is described by the vector 
\ip{t)} = \x; no, ni, . . . ; m) in the Hilbert space H asso- 
ciated with the machine. 

At the moment t = the state of the 
machine is described by the vectors \ip(Q)) = 
£ m a m |0;0,...,0;..., 0,0,0,...) such that 



^ 6(q 1 ,ai,q 2 ,a 2 ,d) = 1. (3) 

(q 2 ,a 2 ,d)eQx Ax {-1,0,1} 



£ki 2 = i- 

i 



(4) 



G 



The evolution of the quantum Turing machine is de- 
scribed by the unitary operator U acting on %. 

A classical probabilistic (or nondeterministic) Turing 
machine can be described as a quantum Turing machine 
such that, at each step of its evolution, the state of the 
machine is represented by the base vector. 

The formal definition of the quantum Turing machine 
was introduced in [21] , 

It is common to use real numbers as amplitudes when 
describing the state of quantum systems during quantum 
computation. To avoid problems with an arbitrary real 
number we introduce the class of numbers which can be 
used as amplitudes for amplitude transition functions of 
the quantum Turing machine. 

Let us denote by C the set of complex numbers c € C, 
such that there exists a deterministic Turing machine, 
which allows to calculate Re (c) and Im (c) with accuracy 
7^- in time polynomial in n. 

Definition 8 (Quantum Turing Machine) A 

quantum Turing machine (QTM) M over an alphabet A 
is a sixtuple (Q,A,S,q Ql q a ,q r ), where 

• Q is the set of internal control states, 

• <Zo i Qa i Qr € Q are initial, accepting and rejecting 
states, 

• 5: QxAxQxAx { — 1, 0, 1} M> C is a transition 
amplitude function i.e. 

^2 \S{qi,a 1 ,q 2 ,a 2 ,d)\ 2 = 1. (5) 

(92,a2,<i)eQx/lx{-l,0a} 

Reversible classical Turing machines (i.e. Turing ma- 
chines with reversible transition function) can be viewed 
as particular examples of quantum machines. Since any 
classical algorithm can be transformed into reversible 
form, it is possible to simulate a classical Turing machine 
using quantum Turing machine. 



4- Quantum complexity 

Quantum Turing machine allows for rigorous analysis 
of algorithms. This is important since the main goal 
of quantum information theory is to provide some gain 
in terms of speed or memory with respect to classical 
algorithms. It should be stressed that at the moment 
no formal proof has been given that a quantum Turing 
machine is more powerful than a classical Turing ma- 
chine [22] . 

In this section we give some results concerning 
quantum complexity theory. See also |21[ 147] for a in- 
troduction to this subject. 

In analogy to classical case it is possible to define com- 
plexity classes for the quantum Turing machine. The 
most important complexity class is this case is BQP. 



Definition 9 Complexity class BQP contains languages 
L for which there exists a quantum Turing machine run- 
ning in polynomial time such that, for any input word x, 
this word is accepted with probability at least | if x G L 
and is rejected with probability at least | if x £ L. 

Class BQP is a quantum counterpart of the classical 
class BPP. 

Definition 10 Complexity class BPP contains lan- 
guages L for which there exists a nondeterministic Turing 
machine running in polynomial time such that, for any 
input word x, this word is accepted with probability at 
least | if x € L and is rejected with probability at least | 
if x $ L. 

Since many results in complexity theory are stated in 
terms of oracles, we define an oracle as follows. 

Definition 11 An oracle or black box is an imaginary 
machine which can decide certain problems in a single 
operation. 

We use notation A B to describe the class of problems 
solvable by an algorithm in class A with an oracle for the 
language B. 

It was shown |21j that the quantum complexity classes 
are related as follows. 

Theorem 2 Complexity classes fulfil the following in- 
equality 

BPP C BQP C P #p . (6) 

Complexity class #P consists of problems of the form 
compute f(x), where f is the number of accepting paths of 
an NP machine. For example problem ^SAT formu- 
lated below is in ffP. 

Problem 1 (#SAT) For a given Boolean formula, 
compute how many satisfying true assignments it has. 

Complexity class P# p consists of all problems solvable 
by a machine running in polynomial time which can use 
oracle for solving problems in #P. 

Complexity ZOO [U] contains the description of com- 
plexity classes and many famous problems from complex- 
ity theory. The complete introduction to the complexity 
theory can be found in [35] . Theory of NP-completeness 
with many examples of problems from this class is presen- 
ted in [48] . 

Many important results and basic definitions concern- 
ing quantum complexity theory can be found in [21] . 
The proof of equivalence between quantum circuit and 
quantum Turing machine was given in [45] . An interest- 
ing discussion of quantum complexity classes and relation 
of BQP class to classical classes can be found in [22'. 
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B. Quantum computational networks 

After presenting the basic facts about Turing machines 
we are ready to introduce more usable models of comput- 
ing devices. We start by defining Boolean circuits and 
extending this model to the quantum case. 

1. Boolean circuits 

Boolean circuits are used to compute functions of the 
form 

/: {0,1}'" ^{0,1}". (7) 

Basic gates (functions) which can be used to define such 
circuits are 

• A : {0,1} 2 i y {0,1}, A(x,y) = l^x = y = l 
(logical and), 

• V : {0,1} 2 i — ^ {0,1}, V(x,y) = 0^x = y = 
(logical or), 

• ~: {0, 1} i v {0, 1}, ~ (x) = 1 — x (logical not). 

The set of gates is called universal if all functions 
{0, 1}™ i-)- {0, 1} can be constructed using the gates from 
this set. It is easy to show that the set of functions com- 
posed of ~, V and A is universal. Thus it is possible 
to compute any functions {0, 1}" i— > {0, l}" 1 using only 
these functions. The full characteristic of universal sets 
of functions was given by Post in 1949 [49] . 

Using the above set of functions a Boolean circuit is 
defined as follows. 

Definition 12 (Boolean circuit) A Boolean circuit is 
an acyclic direct graph with nodes labelled by input vari- 
ables, output variables or logical gates V, A or ~. 

Input variable node has no incoming arrow while output 
variable node has no outcoming arrows. The example of 
a Boolean circuit computing the sum of bits x\ and xi is 
given in Figure [4] 




Figure 4. The example of a Boolean circuit computing the 
sum of bits X\ and X2 [50]. Nodes labelled x\ and X2 represent 
input variables and nodes labelled j/i and yi represent output 
variables. 

Note that in general it is possible to define a Boolean 
circuit using different sets of elementary functions. Since 



functions V , A and ~ provide a universal set of gates we 
defined Boolean circuit using these particular functions. 

Function / : {0, l} m H> {0, 1} is defined on the binary 
string of arbitrary length. Let /„ : {0, l} m i-» {0, 1}™ 
be a restriction of / to {0, 1}™. For each such restriction 
there is a Boolean circuit C n computing /„ . We say that 
Co, C\, C%, . . . is a family of Boolean circuits computing 
/• 

Note that any binary language L C {0, 1}* can be 
accepted by some family of circuits. But since we need 
to know the value of f n to construct a circuit C n such 
family is not an algorithmic device at all. We can state 
that there exists a family accepting the language, but we 
do not know how to build it [35] , 

To show how Boolean circuits are related to Turing 
machines we introduce uniformly generated circuits. 

Definition 13 We say that language L 6 A* has uni- 
formly polynomial circuits if there exists a Turing ma- 
chine M that an input 1 . . . 1 outputs the graph of circuit 

n 

C n using space 0(log?i), and the family Cq,C\,... ac- 
cepts L. 

The following theorem provides a link between uni- 
formly generated circuits and Turing machines. 

Theorem 3 A language L has uniformly polynomial cir- 
cuit iff L e P . 

Quantum circuits model is an analogous to uniformly 
polynomial circuits. They can be introduced as the 
straightforward generalisation of reversible circuits. 

2. Reversible circuits 

The evolution of isolated quantum systems is described 
by a unitary operator U . The main difference with re- 
spect to classical evolution is that this type of evolution 
is reversible. 

Before introducing a quantum circuit we define a re- 
versible Boolean circuit 

Definition 14 (Reversible gate) A classical revers- 
ible function (gate) {0, l} m i— > {0, l} m is a permutation. 

Definition 15 A reversible Boolean circuit is a Boolean 
circuit composed of reversible gates. 

The important fact expressed by the following the- 
orem allows us to simulate any classical computation on 
a quantum machine described using a reversible circuit 

Theorem 4 All Boolean circuits can be simulated using 
reversible Boolean circuits. 

Like in the case of nonreversible circuit one can intro- 
duce the universal set of functions for reversible circuits. 

The important example of a gate universal for revers- 
ible Boolean circuits is a Toffoli gate. The graphical rep- 
resentation of this gate is presented in Figure [5] The 
following theorem was proved by Toffoli [52] . 
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Figure 5. Classical Toffoli gate is universal for reversible cir- 
cuits. It was also used to provide the universal set of quantum 
gates 1511. 



Theorem 5 A Toffoli gate is a universal reversible gate. 

As we will see in the following section it is possible to 
introduce two-bit quantum gates which are universal for 
quantum circuits. This is impossible in the classical case 
and one needs at least a three-bit gate to construct the 
universal set of reversible gates. 

In particular, any reversible circuit is automatically a 
quantum circuit. However, quantum circuits offer much 
more diversity in terms of the number of allowed opera- 
tions. 



3. Quantum circuits 

The computational process of the quantum Turing ma- 
chine is complicated since data as well as control variables 
can be in a superposition of base states. To provide more 
convenient method of describing quantum algorithms one 
can use a quantum circuits model. This model is some- 
times called a quantum gate array model. 

Quantum circuits model was first introduced by 
Deutsch in |51j and it is the most commonly used nota- 
tion for quantum algorithms. It is much easier to ima- 
gine than the quantum Turing machine since the control 
variables (executed steps and their number) are classical. 
There are only quantum data (e.g. qubits or qudits and 
unitary gates) in a quantum circuit. 

A quantum circuit consists of the following elements 
(see Table [n): 

• the finite sequence of wires representing qubits or 
sequences of qubits (quantum registers), 

• quantum gates representing elementary operations 
from the particular set of operations implemented 
on a quantum machine, 

• measurement gates representing a measurement op- 
eration, which is usually executed as the final step 
of a quantum algorithm. It is commonly assumed 
that it is possible to perform the measurement on 
each qubit in canonical basis {|0),|1)} which cor- 
responds to the measurement of the S z observable. 

The concept of a quantum circuit is the natural gener- 
alisation of acyclic logic circuits studied in classical com- 
puter science. Quantum gates have the same number of 



inputs as outputs. Each n-qubit quantum gate repres- 
ents the 2 n -dimensional unitary operation of the group 
SU(2 n ), i.e. generalised rotation in a complex Hilbert 
space. 

The main advantage of this model is its simplicity. It 
also provides very convenient representation of physical 
evolution in quantum systems. 

From the mathematical point of view quantum gates 
are unitary matrices acting on n-dimensional Hilbert 
space. They represent the evolution of an isolated 
quantum system |53j . 

The problem of constructing new quantum algorithms 
requires more careful study of operations used in 
quantum circuit model. In particular wc arc interested in 
efficient decomposition of quantum gates into elementary 
operations. 

We start by providing basic charactersistic of unitary 
matrices l53ll54l 



Theorem 6 Every unitary 2x2 matrix G E U(2) can 
be decomposed using elementary rotations as 



wht 



G = <S>(8)R z (a)R v (6)R z (P) 



<I>,0 = ( B o £ ) 



(8) 



Ry(0 



cos(£/2) sin(£/2)\ 
-sin(£/2) cos(£/2) J ' 



and 



R z (0 



o 



We introduce the definition of quantum gates as stated 
in [50]. 

Definition 16 A quantum gate U acting on m qubits is 
a unitary mapping on C 2 = C 2 <g> . . . <g> C 2 



m times 



u 



(9) 



which operates on the fixed number of qubits. 



Formally, a quantum circuit is defined as the unitary 
mapping which can be decomposed into the sequence of 
elementary gates. 

Definition 17 A quantum circuit on m qubits is a unit- 
ary mapping on C 2 , which can be represented as a con- 
catenation of a finite set of quantum gates. 

Any reversible classical gate is also a quantum gate. 
In particular logical gate ~ (negation) is represented by 
quantum gate NOT, which is realized by a x Pauli matrix. 

As we know any Boolean circuit can be simulated by 
a reversible circuit and thus any function computed by 
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a Boolean circuit can be computed using a quantum cir- 
cuit. Since a quantum circuit operates on a vector in 
complex Hilbcrt space it allows for new operations typ- 
ical for this model. 

The first example of quantum gate which has no clas- 
sical counterpart is \J NOT gate. It has the following 
property 



VNOtVnOT = NOT, 



(10) 



which cannot be fulfilled by any classical Boolean func- 
tion {0,1} i-> {0,1}. Gate \/N is represented by the 
unitary matrix 



Vnot = 



1 



1+i 1 
1-i 1 



(11) 



Another example is Hadamard gate H . This gate is 
used to introduce the superposition of base states. It 
acts on the base state as 

H\0) = 4 (|0) + |1» ,H\1) = ± (|0) - |1» . (12) 



V2 



V2 



If the gate G is a quantum gate acting on one qubit it 
is possible to construct the family of operators acting on 
many qubits. The particularly important class of mul- 
tiqubit operations is the class of controlled operations. 

Definition 18 (Controlled gate) Let G be a 2 x 2 

unitary matrix representing a quantum gate. Operator 



|1)<1|<8G+|0)<0|®I (13) 

acting on two qubits, is called a controlled-G gate. 

Here A (g> B denotes the tensor product of gates (unitary 
operator) A and B, and I is an identity matrix. If in the 
above definition we take G = NOT we get 



|1)<1|®<7*+|0}<0|®I = 



10 
10 
1 
10 



(14) 



which is the definition of CWOT(controlled-Af OT) gate. 
This gate can be used to construct the universal set of 
quantum gates. This gate also allows to introduce en- 
tangled states during computation 



CNOT(H® I) 1 00) 



CNOT^= (|00) 



1 10)) 



V2 



(|00) + |11» 



The classical counterpart of CNOT gate is XOR gate. 

Other examples of single-qubit and two-qubit quantum 
gates are presented in Table [TXJ. In Figure [6] a quantum 
circuit for quantum Fourier transform on three qubits is 
presented. 

One can extend Definition [18] and introduce quantum 
gates with many controlled qubits. 



ko) 
ki> 

l?2> 
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Figure 6. Quantum circuit representing quantum Fourier 
transform for three qubits. Elementary gates used in this 
circuit are described in Table HTl 
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Figure 7. Circuit for quantum teleportation. Double lines 
represent the operation which is executed depending on the 
classical data obtained after the measurement on a subsystem. 



Definition 19 Let G be a 2x2 unitary matrix. Quantum 
gate defined as 

\l^A){l^A\®G+ I0<*l® 1 ( 15 ) 

n-l n-1 1^1... 1 



is called (n - 
by A„_i(G). 



l)-controlled G gate. We denote this gate 



This gate A„_i(G) is sometimes referred to as a gen- 
eralised Toffoli gate or a Toffoli gate with m controlled 
qubits. Graphical representation of this gate is presented 
in Figure [8j 

The important feature of quantum circuits is expressed 
by the following universality property [54] . 

Theorem 7 The set of gates consisting of all one-qubit 
gates U(2) and one two-qubit CNOT gate is universal in 
the sense that any n-qubit operation can be expressed as 
the composition of these gates. 



Xl 


X2 


xi XOR X2 














1 


1 


1 





1 


1 


1 






Table I. Logical values for XOR gate. Quantum CNOT gate 
computes value of x\ XOR X2 in the first register and stores 
values of X2 in the second register. 
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The name of the gate 



Graphical representation 



Mathematical form 



Hadamard 



Pauli X 



Pauli Y 



Pauli Z 



Phase 



H 



1 

^2 



CNOT 



SWAP 

Measurement 
qubit 
n qubits 
classical bit 



X 



Y 



-e- 



i i 

i -i 

oA 

1 ; 

-i \ 

1 / 

1 

-1 

1 ' 
i 



1 

e 17r/i 



/ 1 \ 
10 
1 

V o o i o / 

/ 1 \ 
10 
10 

V o o o i / 
i o \ [oo 

J ' \ 1 
wire = single qubit 



wire representing n qubits 
double wire = single bit 



Table II. Basic gates used in quantum circuits with their graphical representation and mathematical form. Note that measure- 
ment gate is represented in Kraus form, since it is the example of non-unitary quantum evolution. 



Note that, in contrast to the classical case, where one 
needs at least three-bit gates to construct a universal set, 
quantum circuits can be simulated using one two-qubit 
universal gate. 

In order to implement a quantum algorithm one has 
to decompose many qubit quantum gates into element- 
ary gates. It has been shown that almost any n-qubit 
quantum gate (n > 2) can be used to the build a univer- 
sal set of gates |55j in the sense that any unitary opera- 
tion on the arbitrary number of qubits can be expressed 
as the composition of gates from this set. In fact the 
set consisting of two-qubit exclusive-or (XOR) quantum 
gate and all single-qubit gates is also universal [54] . 



Let us assume that we have the set of gates containing 
only CNOT and one-qubit gates. In [SB] theoretical lower 
bound for the number of gates required to simulate a cir- 
cuit using these gates was derived. The efficient method 
of elementary gates sequence synthesis for an arbitrary 
unitary gate was presented in [57] . 

Theorem 8 (Shende-Markov-Bullock) Almost all 
n-qubit operators cannot be simulated by a circuit with 
fewer than \\[A n - 3n - 1]] CNOT gates. 

In [SB] the construction providing the efficient way of 
implementing arbitrary quantum gates was described. 
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Figure 8. Generalised quantum Toffoli gate acting on n 
qubits. Gate G is controlled by the state of n — 1 qubits 
according to Definition |19| 

The resulting circuit has complexity 0(4") which coin- 
cides with lower bound from Theorem [8) 

It is useful to provide more details about the special 
case, when one uses gates with many controlled and one 
target qubits. The following results were proved in [54] . 

Theorem 9 For any single-qubit gate U the gate 
A n _i(£/) can be simulated in terms of 0(n 2 ) basic op- 
erations. 

In many situations it is useful to construct a circuit 
which approximates the required circuit. We say that 
quantum circuits approximate other circuits with ac- 
curacy s if the distance (in terms of Euclidean norm) 
between unitary transformations associated with these 
circuits is at most e [54"] . 

Theorem 10 For any single-qubit gate U and e > 
gate A n -\(U) can be approximated with accuracy e us- 
ing 0(nlog |) basic operations. 

Note that the efficient decomposition of a quantum 
circuit is crucial in physical implementation of quantum 
information processing. In particular case decomposition 
can be optimised using the set of elementary gates spe- 
cific for target architecture. CNOT gates are of big im- 
portance since they allow to introduce entangled states 
during computation. It is also hard to physically realise 
CNOT gate since one needs to control physical interac- 
tion between qubits. 

One should also note that for some classes of quantum 
circuits it is possible to construct their classical counter- 
parts, which can be used to simulate quantum computa- 
tion performed by these circuits efficiently. The most not- 
able class having this property is a class of circuits CHP 
class, which consists of stabilizer circuits, i.e. circuits con- 
sisting solely of CNOT, Hadamard and phase gates [59"] . 
This property is known as so called Gottesman-Knill the- 
orem. 

Theorem 11 (Gottesman-Knill) Any stabilizer cir- 
cuit can be efficiently simulated on a classical machine. 

It is worth noting that gates used to construct sta- 
bilizer circuits do not provide an universal set of gates. 
Nevertheless, such circuits can produce highly entangled 
states. 



C. Random access machines 

Quantum circuit model does not provide a mechan- 
ism for controlling with classical machine the operations 
on quantum memory. Usually quantum algorithms are 
described using mathematical representation, quantum 
circuits and classical algorithms [50]. The model of 
quantum random access machine is built on an assump- 
tion that the quantum computer has to be controlled by 
a classical device [61]. Schematic presentation of such 
architecture is provided in Figure [9] 

Quantum random access machine is interesting for 
us since it provides a convenient model for developing 
quantum programming languages. However, these lan- 
guages are our main area of interest. We see no point 
in providing the detailed description of this model as it 
is given in [61] together with the description of hybrid 
architecture used in quantum programming. 



1. Classical RAM model 

The classical model of random access machine (RAM) 
is the example of more general register machines [361 137] . 

The random access machine consists of an unboun- 
ded sequence of memory registers and a finite number of 
arithmetic registers. Each register may hold an arbitrary 
integer number. The programme for the RAM is a finite 
sequence of instructions II = (tti, . . . , n n ). At each step 
of execution register i holds an integer and the ma- 
chine executes instruction 7r K , where n is the value of the 
programme counter. Arithmetic operations are allowed 
to compute the address of a memory register. 

Despite the difference in the construction between a 
Turing machine and RAM, it can be easily shown that 
a Turing machine can simulate any RAM machine with 
polynomial slow-down only [35] . 

It is worth noting that programming languages can 
be defined without using RAM model. Interesting pro- 
gramming language for a Turing machine V" , providing 
the minimal set of instructions, was introduced by Bohm 
in [43] . 



2. Quantum RAM model 

Quantum random access machine (QRAM) model is 
the extension of the classical RAM. QRAM can exploit 
quantum resources and, at the same time, can be used 
to perform any kind of classical computation. It allows 
us to control operations performed on quantum registers 
and provides the set of instructions for defining them. 

Recently a new model of sequential quantum random 
machine (SQRAM) has been proposed. Instruction set 
for this model and compilation of high-level languages is 
discussed in [52]. However, it is very similar to QRAM 
model. 
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the sequence of elementary gates 

the description of 




Figure 9. The model of classically controlled quantum ma- 
chine (61J . Classical computer is responsible for perform- 
ing unitary operations on quantum memory. The results of 
quantum computation are received in the form of measure- 
ment results. 

The quantum part of QRAM model is used to generate 
probability distribution. This is achieved by performing 
measurement on quantum registers. The obtained prob- 
ability distribution has to be analysed using a classical 
computer. 



Procedure: FouRlER(a, d) 

Input: A quantum register a with d qubits. Qubits are- 
numbered from to d — 1. 

Output: The amplitudes of a are Fourier transformed 
over Z 2 <j . 

C: assign value to classical variable 
u *- e l2 ^ 2d 

C: perform sequence of gates 
for i = d— 1 to i = 
for j = d — 1 to j — i + 1 
if aj then lZ^ 2 d-i-i+j (ch) 
C: number of loops executing phase 
C: depends on the required accuracy 
C: of the procedure 
H{ai) 

C: change the order of qubits 
for j = to j = i - 1 
SWAP(aj, a d -a-j ) 

Listing 1. Quantum pseudocode for quantum Fourier trans- 
form on d qubits. Quantum circuit for this operation with 
d = 3 is presented in Figure [6] 



3. Quantum pseudocode 

Quantum algorithms are, in most of the cases, de- 
scribed using the mixture of quantum gates, mathemat- 
ical formulas and classical algorithms. The first attempt 
to provide a uniform method of describing quantum al- 
gorithms was made in |63| . where the author introduced 
a high-level notation based on the notation known from 
computer science textbooks [64] . 

In [BU] the first formalised language for description of 
quantum algorithms was introduced. Moreover, it was 
tightly connected with the model of quantum machine 
called quantum random access machine (QRAM). 

Quantum pseudocode proposed by Knill [60] is based 
on conventions for classical pseudocode proposed in [Ml 
Chapter 1]. Classical pseudocode was designed to be 
readable by professional programmers, as well as people 
who had done a little programming. Quantum pseudo- 
code introduces operations on quantum registers. It also 
allows to distinguish between classical and quantum re- 
gisters. 

Quantum registers are distinguished by underlining 
them. They can be introduced by applying quantum op- 
erations to classical registers or by calling a subroutine 
which returns a quantum state. In order to convert a 
quantum register into a classical register measurement 
operation has to be perfomed. 

The example of quantum pseudocode is presented in 
Listing [l] It shows the main advantage of QRAM model 
over quantum circuits model - the ability to incorpor- 
ate classical control into the description of quantum al- 
gorithm. 

Operation rl(ai) executes a quantum Hadamard gate 



on a quantum register a; and SWAV(a,i, aj) performs 
SWAP gate between and aj. Operation lZ$(ai) that 
executes a quantum gate R((f>) is defined as 



R{4>) 




(16) 



on the quantum register a^. Using conditional contrac- 
tion 

if aj then TZ^ica) 

it is easy to define controlled phase shift gate (see Defin- 
19 1. Similar construction exists in QCL quantum 



it ion 



programming language described in Section III 



The measurement of a quantum register can be indic- 
ated using an assignement a, aj. 



4- Quantum programming environment 

Since the main aim of this paper is to present the ad- 
vantages and limitations of high-level quantum program- 
ming languages, we need to explain how these languages 
are related to quantum random access machine. Thus 
as the summary of this section we present the overview 
of an architecture for quantum programming, which is 
based on the QRAM model. 

The architecture proposed in [BSl [SB] is designed for 
transforming a high-level quantum programming lan- 
guage into the technology-specific implementation set of 
operations. This architecture is composed of four layers: 

• High level programming language provid- 
ing high-level mechanisms for performing useful 



13 



quantum computation; this language should be in- 
dependent from particular physical implementation 
of quantum computing. 

• Compiler of this language providing architec- 
ture independent optimisation; also compilation 
phase can be used to handle quantum error cor- 
rection required to perform useful quantum com- 
putation. 

• Quantum assembly language (QASM) - as- 
sembly language extended by the set of instructions 
used in the quantum circuit model. 

• Quantum physical operations language (QC- 
POL), which describes the execution of quantum 
programme in a hardware-dependent way; it in- 
cludes physical operations and it operates on the 
universal set of gates optimal for a given physical 
implementation. 

The authors of [SS] do not define a specific high- 
level quantum programming language. They point out, 
however, that existing languages, mostly based on Dirac 
notation, do not provide the sufficient level of abstrac- 
tion. They also stress, following [ST], that it should 
have the basic set of features. We will discuss these 
basic requirements in detail in Section |III| At the mo- 
ment quantum assembly language (QASM) is the most 
interesting part of this architecture, since it is tightly 
connected to the QRAM model. 

QASM should be powerful enough for representing 
high level quantum programming language and it should 
allow for describing any quantum circuit. At the same 
time it must be implementation-independent so that it 
could be used to optimise the execution of the programme 
with respect to different architectures. 

QASM uses qubits and cbits (classical bit) as basic 
units of information. Quantum operations consist of 
unitary operations and measurement. Moreover, each 
unitary operator is expressed in terms of single qubit 
gates and CNOT gates. 

In the architecture proposed in [53] each single-qubit 
operation is stored as the triple of rationals. Each ra- 
tional multiplied by 7r represents one of three Eulcr- 
angles, which are sufficient to specify one-qubit opera- 
tion. 



III. QUANTUM PROGRAMMING 
LANGUAGES 

Quantum algorithms [12 HU EH1 EH! and communica- 
tion protocols [6j [70l [71] are described using a language 
of quantum circuits |53j . While this method is conveni- 
ent in the case of simple algorithms, it is very hard to 
operate on compound or abstract data types like arrays 
or integers using this notation |T9l [72] . 

This lack of data types and control structures motiv- 
ated the development of quantum pseudocode [BUI [73] 



and various quantum programming languages 61, 66, 74- 

EZj. 

Several languages and formal models were proposed for 
the description of quantum computation process. The 
most popular of them is quantum circuit model [51] . 
which is tightly connected to the physical operations im- 
plemented in the laboratory. On the other hand the 
model of quantum Turing machine is used for analysing 
the complexity of quantum algorithms [21] . 

Another model used to describe quantum computers 
is Quantum Random Access Machine (QRAM). In this 
model we have strictly distinguished the quantum part 
performing computation and the classical part, which is 
used to control computation. This model is used as a 
basis for most quantum programming languages [781480] . 
Among high-level programming languages designed for 
quantum computers we can distinguish imperative and 
functional languages. 

At the moment of writing this paper the most advanced 
imperative quantum programming language is Quantum 
Computation Language (QCL) designed and implemen- 
ted by Omer [5TJ [STJ [SS] . QCL is based on the syntax 
of C programming language and provides many elements 
known from classical programming languages. The in- 
terpreter is implemented using simulation library for ex- 
ecuting quantum programmes on classical computer, but 
it can be in principle used as a code generator for classical 
machine controlling a quantum circuit. 

Along with QCL several other imperative quantum 
programming languages were proposed. Notably Q Lan- 
guage developed by Betteli [67j [74] and libquantum [83] 
have the ability to simulate noisy environment. Thus, 
they can be used to study decoherence and analyse the 
impact of imperfections in quantum systems on the ac- 
curacy of quantum algorithms. 

Q Language [84] is implemented as a class library for 
C++ programming language and libquantum is imple- 
mented as a C programming language library. Q Lan- 
guage provides classes for basic quantum operations like 
QHadamard, QFourier, QNot, QSwap, which are derived 
from the base class Qop. New operators can be defined 
using C++ class mechanism. Both Q Language and 
libquantum share some limitation with QCL, since it is 
possible to operate on single qubits or quantum registers 
(i.e. arrays of qubits) only. Thus, they are similar to 
packages for computer algebra systems used to simulate 
quantum computation |85[ 186] . 

Concerning problems with physical implementations 
of quantum computers, it became clear that one needs 
to take quantum errors into account when modelling 
quantum computational process. Also quantum com- 
munication has become very promising application of 
quantum information theory over the last few years. 
Both facts are reflected in the design of new quantum 
programming languages. 

LanQ developed by Mlnaffk was defined in [77l [87]. 
It provides syntax based on C programming language. 
LanQ provides several mechanisms such as the creation 
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of a new process by forking and interprocess communic- 
ation, which support the implementation of multi-party 
protocols. Moreover, operational semantics of LanQ has 
been defined. Thus, it can be used for the formal reas- 
oning about quantum algorithms. 

It is also worth to mention new quantum program- 
ming languages based on functional paradigm. Research 
in functional quantum programming languages started 
by introducing quantum lambda calculus f88 . It was in- 
troduced in a form of simulation library for Scheme pro- 
gramming language. QPL [89] was the first functional 
quantum programming language. This language is stat- 
ically typed and allows to detect errors at compile-time 
rather than run-time. 

A more mature version of QPL is cQPL — communic- 
ation capable QPL [75]. cQPL was created to facilitate 
the development of new quantum communication proto- 
cols. Its interpreter uses QCL as a backend language so 
cQPL programmes are translated into C++ code using 
QCL simulation library. 



A. Requirements for quantum programming 
language 

Taking into account QRAM model described in Section 



Table III contains the comparison of several quantum 
programming languages. It includes the most import- 
ant features of existing languages. In particular we list 
the underlying mathematical model (i.e. pure or mixed 
states) and the support for quantum communication. 



All languages listed in Table III are universal and thus 
they can be used to compute any function computable on 
a quantum Turing machine. Consequently all these lan- 
guage provide the model of quantum computation which 
is equivalent to the model of a quantum Turing machine. 

In this section we compare the selected quantum 
programming languages and provide some examples of 
quantum algorithms and protocols implemented in these 
languages. We also describe their main advantages and 
limitations. We introduce the basic syntax of three of the 
languages listed in Table |III| - QCL, LanQ and cQPL. 
This is motivated by the fact that these languages have 
a working interpreter and can be used to perform simu- 
lations of quantum algorithms. We introduce basic ele- 
ments of QCL required to understand basic programmes. 
We also compare the main features of the presented lan- 
guages. 

The main problem with current quantum programming 
languages is that they tend to operate on very low-level 
structures only. In QCL quantum memory can be ac- 
cessed using only qreg data type, which represents the 
array of qubits. In the syntax of cQPL data type qint 
has been introduced, but it is only synonymous for the 
array of 16 qubits. Similar situation exists in LanQ |87j . 
where quantum data types are introduced using qnit 
keyword, where n represents a dimension of elementary 
unit (e.g. for qubits n — 2, for qutrits n — 3). How- 
ever, only unitary evolution and measurement can be 
performed on variables defined using one of these types. 



II C 2 we can formulate basic requirements which have to 
be fulfilled by any quantum programming language |74j . 

• Completeness: Language must allow to express 
any quantum circuit and thus enable the program- 
mer to code every valid quantum programme writ- 
ten as a quantum circuit. 

• Extensibility: Language must include, as its sub- 
set, the language implementing some high level 
classical computing paradigm. This is important 
since some parts of quantum algorithms (for ex- 
ample Shor's algorithm) require nontrivial classical 
computation. 

• Separability: Quantum and classical parts of the 
language should be separated. This allows to ex- 
ecute any classical computation on purely classical 
machine without using any quantum resources. 

• Expressivity: Language has to provide high level 
elements for facilitating the quantum algorithms 
coding. 

• Independence: The language must be independ- 
ent from any particular physical implementation of 
a quantum machine. It should be possible to com- 
pile a given programme for different architectures 
without introducing any changes in its source code. 

As we will see, the languages presented in this Section 
fulfil most of the above requirements. The main problem 
is the expressivity requirement. 



B. Imperative quantum programming 

First we focus on quantum programming languages 
which are based on the imperative paradigm. They in- 
clude quantum pseudocode, discussed in Section II C 2[ 
Quantum Computation Language (QCL) created by 
Omer [5TJ [S3 [52] and LanQ developed by Mlnank 

[7J [771 M\ 

Below we provide an introduction to QCL. It is one 
of the most popular quantum programming languages. 
Next, we introduce the basic elements of LanQ. This lan- 
guage provides the support for quantum protocols. This 
fact reflects the recent progress in quantum communica- 
tion theory. 



1. Quantum Computation Language 

QCL (Quantum Computation Language) [5TJ [5TJ [52] is 
the most advanced implemented quantum programming 
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Table III. The comparison of quantum programming languages with information about implementation and basic features. 
Based on information from [75] and |87| . 



language. Its syntax resembles the syntax of C program- 
ming language ^90] and classical data types are similar to 
data types in C or Pascal. 

The basic built-in quantum data type in QCL is qureg 
(quantum register) . It can be interpreted as the array of 
qubits (quantum bits). 

qureg xl[2]; // 2—qubit quantum register xl 
qureg x2[2]; // 2—qubit quantum register x2 
H(xl); // Hadamard operation on xl 
H(x2[l]); // and on the second qubit of x2 

Listing 2. Basic operations on quantum registers and subre- 
gisters in QCL. 

QCL standard library provides standard quantum op- 
erators used in quantum algorithms, such as: 

• Hadamard H and Not operations on many qubits, 

• controlled not (CNot) with many target qubits and 
Swap gate, 

• rotations: RotX, RotY and RotZ, 

• phase (Phase) and controlled phase (CPhase). 

Most of them are described in Table Hi] in Section III C 21 
Since QCL interpreter uses qlib simulation library, it 
is possible to observe the internal state of the quantum 
machine during the execution of quantum programmes. 
The following sequence of commands defines two-qubit 
registers a and b and executes H and CNot gates on these 
registers. 

qcl> qureg a [2] ; 
qcl> qureg b [2] ; 
qcl> H(a) ; 

[4/32] 0.5 |0,0> + 0.5 |1,0> + 

0.5 |2,0> + 0.5 |3,0> 
qcl> dump 

: STATE: 4/32 qubits allocated , 
28 / 32 qubits free 

0.5 |0> + 0.5 |1> + 0.5 |2> + 0.5 |3> 
qcl> CNot(a [1] ,b) 

[4/32] 0.5 1 ,0 > + 0.5 |1,0> + 0.5 |2,0> 
+ 0.5 1 3 ,0 > 
qcl> dump 

: STATE: 4/32 qubits allocated , 
28 / 32 qubits free 

0.5 |0> + 0.5 |1> + 0.5 |2> + 0.5 |3> 



Using dump command it is possible to inspect the in- 
ternal state of a quantum computer. This can be help- 
ful for checking if our algorithm changes the state of 
quantum computer in the requested way. 

One should note that dump operation is different from 
measurement, since it does not influence the state of 
quantum machine. This operation can be realised using 
simulator only. 

a. Quantum memory management Quantum 
memory can be controlled using quantum types qureg, 
quconst, quvoid and quscratch. Type qureg is used as 
a base type for general quantum registers. Other types 
allow for the optimisation of generated quantum circuit. 
The summary of types defined in QCL is presented in 
Table Jv] 

b. Classical and quantum procedures and functions 
QCL supports user-defined operators and functions 
known from languages like C or Pascal. Classical sub- 
routines are defined using procedure keyword. Also 
standard elements, known from C programming lan- 
guage, like looping (e.g. for i=l ton] ... }) and con- 
ditional structures (e.g. if x==0 { ... }), can be used to 
control the execution of quantum and classical elements. 
In addition to this, it provides two types of quantum 
subroutines. 

The first type is used for unitary operators. Using 
it one can define new operations, which in turn can be 
used to manipulate quantum data. For example operator 
diffuse defined in Listing[3] defines inverse about the mean 
operator used in Grover's algorithm [14]. This allows to 
define algorithms on the higher level of abstraction and 
extend the library of functions available for a program- 
mer. 

Using subroutines it is easy to describe quantum al- 
gorithms. Figure [10] presents QCL implementation of 
Deutsch's algorithm, along with the quantum circuit for 
this algorithm. This simple algorithm uses all main ele- 
ments of QCL. It also illustrates all main ingredients of 
existing quantum algorithms. 

The second type of quantum subroutine is called a 
quantum function. Quantum functions are also called 
pseudo-classic operators. It can be defined using qufunct 
keyword. The subroutine of type qufunct is used for all 
transformations of the form 

|n> = |/(n)), (17) 
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Type 


Description 


Usage 


qureg 
quvoid 
quconst 
quscratch 


general quantum register 

register which has to be empty when operator is called 

must be invariant for all operators used in quantum conditions 

register which has to be empty before and after the operator is 
called 


basic type 
target register 
quantum conditions 
temporary registers 



Table IV. Types of quantum registers used for memory management in QCL. 



where |n) is a base state and / is a one-to-one Boolean 
function. The example of quantum function is presented 
in Listing [5] 

c. Quantum conditions QCL introduces quantum 
conditional statements, i.e. conditional constructions 
where quantum state can be used as a condition. 

QCL, as well as many classical programming lan- 
guages, provides the conditional construction of the form 

if be then 

block 

where be is a Boolean expression and block is a sequence 
of statements. 

QCL provides the means for using quantum variables 
as conditions. Instead of a classical Boolean variable, the 
variable used in condition can be a quantum register. 

qureg a [2] ; 
qureg b [ 2 ] ; 

// the sequence of statements 
// ... 

// perform CNot if a=|l...l) 
if a { 

CNot(b[0] , b [1] ) ; 

} 

Listing 4. Example of a quantum conditional statement in 
QCL 

In this situation QCL interpreter builds and executes 
the sequence of CNOT gates equivalent to the above 
condition. Here register a is called enable register. 

In addition, quantum conditional structures can be 
used in quantum subroutines. Quantum operators and 
functions can be declared as conditional using cond 
keyword. For example 

operator diffuse (qureg q) { 

H(q); // Hadamard Transform 

Not(q); // Invert q 

CPhase ( pi , q ) ; // Rotate if q = llll.. 

!Not(q); // undo inversion 

!H(q); // undo Hadamard Transform 

} 

Listing 3. The implementation of the inverse about the mean 
operation in QCL [61]. Constant pi represents number 7r. 
Exclamation mark ! is used to indicate that the interpreter 
should use the inverse of a given operator. Operation diffuse 
is used in the quantum search algorithm |14j . 



// conditional phase gate 

extern cond operator Phase (real phi); 

// conditional not gate 

extern cond qufunct Not (qureg q); 



\x) — H- 







H 








H 


\y) 







operator U(qureg x . qureg y) { 

H(x); 

Oracle (x , y ) ; 
H(x & y); 

} 

// classical control structure 
procedure deutsch() { 

// allocate 2 qubits 

qureg x [ 1 ] ; 

qureg y [ 1 ] ; 

int m; 

// evaluation loop 
{ 

// initialise machine state 
reset ; 

// do unitary computation 
U(x,y); 

// measure 2nd register 
measure y ,m; 

// value in 1st register valid? 
} until m==l; 

// measure 1st register which 
measure x ,m; 

// contains g(0) xor g(l) 
print "g(0) „xor„g(l) u =" ,m; 
// clean up 
reset ; 

} 

Figure 10. Quantum circuit for Deutsch's algorithm and QCL 
implementation of this algorithm (see |61| for more examples) . 
Evaluation loop is composed of preparation (performed by 
reset instruction), unitary evolution (U(x,y) operator) and 
measurement. Subroutine Oracle() implements function used 
in Deutsch's algorithm [4"21 191] . 
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declares a conditional Phase gate and a controlled NOT 
gate. Keyword extern indicates that the definition of 
a subroutine is specified in an external file. The enable 
register (i.e. quantum condition) is passed as an implicit 
parameter if the operator is used within the body of a 
quantum if- statement. 

// increment register 
cond qufunct inc(qureg x) { 
int i ; 

for i = #x— 1 to step —1 { 

// apply controlled — not from MSB to LSB 
CNot(x[i] ,x[0:: i ]); 

} 

} 

// equivalent implementation 
// with constant enable register 
// conditional increment as selection operator 
qufunct cinc(qureg x , quconst e) { 
int i ; 

for i = #x— 1 to step —1 { 

CNot(x[i] ,x[0::i]&e); 

} 

} 

Listing 5. Operator for incrementing quantum state in QCL 
defined as a conditional quantum function. Subroutine inc is 
defined using cond keyword and does not require the second 
argument of type quconst. Subroutine cine provides equi- 
valent implementation with explicit-declared enable register. 

In the case of inc procedure, presented in Listing[5j the 
enable register is passed as an implicit argument. This 
argument is set by a quantum if-statement and trans- 
parently passed on to all suboperators. As a result, all 
suboperators have to be conditional. This is illustrated 
by the following example [61j 

// counting and control registers 
qcl> qureg q[4];qureg e[l]; 
// prepare test state 
qcl> H(q[3] & e); 

[5/32] 0.5 |0,0> + 0.5 |8,0> + 0.5 |0,1> 
+ 0.5 |8,1> 

// conditional increment 
qcl> cine (q , e ) ; 

[5/32] 0.5 1 ,0 > + 0.5 |8,0> + 0.5 |1,1> 
+ 0.5 |9,1> 

// equivalent to cinc(q,e) 
qcl> if e { inc(q); } 

[5/32] 0.5 |0,0> + 0.5 |8,0> + 0.5 |2,1> 
+ 0.5 | 10 ,1 > 

// conditional decrement 
qcl> ! cine (q , e ) ; 

[5/32] 0.5 |0,0> + 0.5 |8,0> + 0.5 |1,1> 
+ 0.5 |9,1> 

// equivalent to ! cine ( q , e ) ; 
qcl> if e { ! inc (q) ; } 

[5/32] 0.5 j0,0> + 0.5 |8,0> + 0.5 |0,1> 
+ 0.5 |8,1> 

Finally we should note that a conditional subroutine 
can be called outside a quantum if-statement. In such 



situation enable register is empty and, as such, ignored. 
Subroutine call is in this case unconditional. 



2. LanQ 

Imperative language LanQ is the first quantum pro- 
gramming language with full operation semantics spe- 
cified [87]. 

Its main feature is the support for creating multipart- 
ite quantum protocols. LanQ, as well as cQPL presented 
in the next section, are built with quantum communic- 
ation in mind. Thus, in contrast to QCL, they provide 
the features for facilitating simulation of quantum com- 
munication. 

Syntax of the LanQ programming language is very sim- 
ilar to the syntax of C programming language. In par- 
ticular it supports: 

• Classical data types: int and void. 

• Conditional statements of the form 
if ( cond ) { 

} else { 
} ' 

• Looping with while keyword 
while ( cond ) { 

} ' 

• User-defined functions, for example 

int fun ( int i ) { 
int res ; 

return res ; 

} 

a. Process creation LanQ is built around the con- 
cepts of process and interprocess communication, known 
for example from UNIX operating system. It provides the 
support for controlling quantum communication between 
many parties. The implementation of teleportation pro- 
tocol presented in Listing [6] provides an example of LanQ 
features, which can be used to describe quantum commu- 
nication. 

Function mainQ in Listing [6] is responssible for con- 
trolling quantum computation. The execution of pro- 
tocol is divided into the following steps: 

1. Creation of the classical channel for communicating 
the results of measurement: 

channelpnt] c withends [c0,cl];. 

2. Creation of Bell state used as a quantum channel 
for teleporting a quantum state 

(psiEPR aliasfor [psil, psi2]); this is accomplished 
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by calling external function createEPR() creating an 
entangled state. 

3. Instruction fork executes alice () function, which is 
used to implement sender; original process contin- 
ues to run. 

4. In the last step function bob() implementing a re- 
ceiver is called. 



void al i c e (channelEnd [ int ] cO , 

qbit auxTeleportState ) { 

int i ; 
qbit phi ; 

// prepare state to be teleported 

phi = computeSomething ( ) ; 

// Bell measurement 

i = measure (BellBasis , phi , 

auxTeleportState ) ; 

send ( cO , i ) : 

} 

void bob (channelEnd [ int ] cl , 

qbit stateToTeleportOn ) { 

int i ; 

i = recv ( cl ) ; 

// execute one of the Pauli gates 
// according to the protocol 
if (i = 1) { 

Sigma_z (stateToTeleportOn ) ; 
} else if (i = 2) { 

Sigma_x ( stateToTeleportOn ) ; 
} else if (i = 3) { 

Sigma_x ( stateToTeleportOn ) ; 

Sigma_z( stateToTeleportOn); 

} 

dump_q( stateToTeleportOn ) ; 

} 

void main() { 

channel [ int ] c withends [cO,cl]; 

qbit psil , psi2 ; 

psiEPR aliasfor [psil , psi2 ] ; 

psiEPR = createEPRQ; 

c = new channel [ int ] ( ) ; 

fork alice(cO, psil); 
bob ( cl , psi2 ) ; 

} 

Listing 6. Teleportation protocol implemented in LanQ [87] . 
Functions Sigma_x(), Sigma_y() and Sigma_z() are respons- 
ible for implementing Pauli matrices. Function createEPR() 
(not defined in the listing) creates maximally entangled state 
between parties — Alice and Bob. Quantum communication 
is possible by using the state, which is stored in a global vari- 
able psiEPR. Function computeSomething() (not defined in 
the listing) is responsible for preparing a state to be teleported 
by Alice. 



b. Communication Communication between parties 
is supported by providing send and recv keywords. Com- 
munication is synchronous, i.e. recv delays programme 
execution until there is a value received from the channel 
and send delays a programme run until the sent value is 
received. 

Processes can allocate channels. It should be stressed 
that the notion of channels used in quantum program- 
ming is different from the one used in quantum mech- 
anics. In quantum programming a channel refers to a 
variable shared between processes. In quantum mechan- 
ics a channel refers to any quantum operation. 

Another feature used in quantum communication is 
variable aliasing. In the teleportation protocol presented 
in Listing [6] the syntax for variable aliasing 

qbit psil , psi2 ; 

psiEPR aliasfor [psil , psi2 ] ; 

is used to create quantum state shared among two 
parties. 

c. Types Types in LanQ are used to control the sep- 
aration between classical and quantum computation. In 
particular they are used to prohibit copying of quantum 
registers. The language distinguishes two groups of vari- 
ables [EH Chapter 5]: 

• Duplicable or non-linear types for representing clas- 
sical values, e.g. bit, int, boolean. The value of a 
duplicable type can be exactly copied. 

• Non-duplicable or linear types for controlling 
quantum memory and quantum resources, e.g. 
qbit, qtrit channels and channel ends (see example 



in Listing \6 ) . 
for cloning 



Types from this group do not allow 



One should note that quantum types defined in LanQ 
are mainly used to check validity of the program before its 
run. However, such types do not help to define abstract 
operations. As a result, even simple arithmetic opera- 
tions have to be implemented using elementary quantum 
gates, e.g. using quantum circuits introduced in |93| . 



C. Functional quantum programming — QPL and 
cQPL 

During the last few years few quantum programming 
languages based on functional programming paradigm 
have been proposed [94]. As we have already point out, 
the lack of progress in creating new quantum algorithms 
is caused by the problems with operating on complex 
quantum states. Classical functional programming lan- 
guages have many features which allow to clearly express 
algorithms 41j. In particular they allow for writing bet- 
ter modularised programmes than in the case of imperat- 
ive programming languages 95 . This is important since 
this allows to debug programmes more easily and re- 
use software components, especially in large and complex 
software projects. 
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Quantum functional programming attempts to merge 
the concepts known from classical function programming 
with quantum mechanics. The program in functional 
programming language is written as a function, which 
is defined in terms of other functions. Classical func- 
tional programming languages contain no assignment 
statements, and this allows to eliminate side-effects. [9^] 
It means that function call can have no effect other than 
to compute its result [95] . In particular it cannot change 
the value of a global variable. 

The first attempts to define a functional quantum pro- 
gramming language were made by using quantum lambda 
calculus [85], which was based on lambda calculus. For 
the sake of completeness we can also point out some re- 
search on modelling quantum computation using Haskell 
programming language [971 198] . However, here we focus 
on high-level quantum programming languages. Below 
we present recently proposed languages QPL and cQPL, 
which are based on functional paradigm. They aim to 
provide mechanisms known from programming languages 
like Haskell [53] to facilitate the modelling of quantum 
computation and quantum communication. 

In [Sn] Quantum Programming Language (QPL) was 
described and in [75] its extension useful for modelling 
of quantum communication was proposed. This exten- 
ded language was name cQPL - communication capable 
QPL. Since cQPL compiler is also QPL compiler, we will 
describe cQPL only. 

The compiler for cQPL language described in [75] is 
built on the top of libqc simulation library used in QCL 
interpreter. As a result, cQPL provides some features 
known from QCL. 

Classical elements of cQPL are very similar to classical 
elements of QCL and LanQ. In particular cQPL provides 
conditional structures and loops introduced with while 
keyword. 

new int loop := 10; 
while (loop > 5) do { 

print loop ; 

loop := loop — 1; 

}; 

if (loop = 3) then { 

print " loop „ i s „equal „3" ; 
} else { 

print " loop u i s „not „equal „3" ; 

}; 

Listing 7. Classical control structures in cQPL. 



Procedures 

Procedures can be defined to improve modularity of 
programmes. 

proc test: a : int , q:qbit { 
} ' 



Procedure call has to know the number of parameters re- 
turned by the procedure. If, for example, procedure test 
is defined as above, it is possible to gather the calculated 
results 

new int al = 0; 
new int cv = ; 
new int qv = ; 

(al) := call test(cv, qv); 

or ignore them 

call t est (cv , qv); 

In the first case the procedure returns the values of input 
variables calculated at the end of its execution. 

Classical variables are passed by value i.e. their value 
is copied. This is impossible for quantum variable, since 
a quantum state cannot be cloned [92]. Thus, it is also 
impossible to assign the value of quantum variable cal- 
culated by procedure. 

Note that no cloning theorem requires quantum vari- 
ables to be global. This shows that in quantum case it 
is impossible to avoid some effects known from imperat- 
ive programming and typically not present in functional 
programming languages. 

Global quantum variables are used in Listing [9] to 
create a maximally entangled state in a teleportation 
protocol. Procedure createEPR(eprl, epr2) operates on 
two quantum variables (subsystems) and produces a Bell 
state. 



Quantum elements 

Quantum memory can be accessed in cQPL using vari- 
ables of type qbit or qint. Basic operations on quantum 
registers are presented in Listing [8j In particular, the 
execution of quantum gates is performed by using *= 
operator. 

new qbit ql := 0; 
new qbit q2 := 1; 

// execute CNOT gate on both qubits 
ql , q2 *= CNot ; 

// execute phase gate on the first qubit 
ql *= Phase 0.5; 

Listing 8. State initialisation and basic gates in cQPL. Data 
type qbit represents a single qubit. 

It should be pointed out that qint data type provides 
only a shortcut for accessing the table of qubits. 

Only a few elementary quantum gates are built into 
the language: 

• Single qubit gates H, Phase and NOT implementing 
basic gates listed in Table |H| in Section [II C 2| 

• CNDT operator implementing controlled negation 
and FT(n) operator for n-qubit quantum Fourier 
transform. 
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This allows to simulate an arbitrary quantum computa- 
tion. Besides, it is possible to define gates by specifying 
their matrix elements. 

Measurement is performed using measure/then 
keywords and print command allows to display the value 
of a variable. 

measure a then { 

print "a„is „|0>" ; 
} else { 

print "a„is „|1>" ; 

}; 

In similar manner like in QCL, it is also possible to 
inspect the value of a state vector using dump command. 



Quantum communication 

The main feature of cQPL is its ability to build and 
test quantum communication protocols easily. Commu- 
nicating parties are described using modules. In analogy 
to LanQ, cQPL introduces channels, which can be used 
to send quantum data. Once again we stress that notion 
of channels used in cQPL and LanQ is different from 
that used in quantum theory. Quantum mechanics intro- 
duces channels to describe allowed physical transforma- 
tions, while in quantum programming they are used to 
describe communication links. 

Communicating parties are described by modules, in- 
troduced using module keyword. Modules can exchange 
quantum data (states). This process is accomplished us- 
ing send and receive keywords. 

To compare cQPL and LanQ one can use the imple- 
mentation of the teleportation protocol. The implement- 
ation of teleportation protocol in cQPL is presented in 
Listing [9] while the implementation in LanQ is provided 
in Listing [6] 

module Alice { 

proc createEPR : a:qbit, b:qbit { 

a *= H; 

b , a *= CNot ; 

/* b : Control, a: Target */ 
} in { 



measurement results to Bob */ 
send ml , m2 to Bob ; 



}; 



new qbit 
new qbit 
new qbit 



t e lep ort 
eprl := 
epr2 := 



0; 



call createEPR ( eprl , epr2 ) ; 
send epr2 to Bob; 

/* teleport: Control , eprl: Target 

(see: Figure [?) */ 
teleport , eprl *= CNot; 

new bit ml := 0; 

new bit m2 := 0; 

ml := measure teleport ; 

m2 := measure eprl ; 

/* Transmit the classical 



module Bob { 

receive q:qbit from Alice; 
receive ml: bit, m2:bit from Bob; 

if (ml = 1) then { q *= [[ 0,1,1,0 ]]; 
/* Apply sigma_x */ }; 

if (m2 = 1) then { q *= [[ 1,0,0,-1 ]] ; 
/* Apply sigma_z */}; 

/* The state is now telep orted */ 
dump q ; 

}; 

Listing 9. Teleportation protocol implemented in cQPL (from 
[75]). Two parties - Alice and Bob - are described by 
modules. Modules in cQPL are introduced using module 
keyword. 



IV. SUMMARY 

The main goal of this paper is to acquaint the reader 
with quantum programming languages and computa- 
tional models used in quantum information theory. We 
have described a quantum Turing machine, quantum cir- 
cuits and QRAM models of quantum computation. We 
have also presented three quantum programming lan- 
guages - namely QCL, LanQ and cQPL. 

First we should note that the languages presented in 
this paper provide very similar set of basic quantum gates 
and allow to operate only on the arrays of qubits. Most of 
the gates provided by these languages correspond to the 
basic quantum gates presented in Section |II C 2| Thus, 
one can conclude that the presented languages have the 
ability to express quantum algorithms similar to the abil- 
ities of a quantum circuit model. 

The biggest advantage of quantum programming lan- 
guages is their ability to use classical control structures 
for controlling the execution of quantum operators. This 
is hard to achieve in quantum circuits model and it re- 
quires the introduction of non-unitary operations to this 
model. In addition, LanQ and cQPL provide the syntax 
for clear description of communication protocols. 

The syntax of presented languages resembles the syn- 
tax of popular classical programming languages from the 
C programming language family |90) . As such, it can 
be easily mastered by programmers familiar with clas- 
sical languages. Moreover, the description of quantum 
algorithms in quantum programming languages is bet- 
ter suited for people unfamiliar with the notion used in 
quantum mechanics. 

The main disadvantage of described languages is the 
lack of quantum data types. The types defined in de- 
scribed languages are used mainly for two purposes: 
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• To avoid compile-time errors caused by copying of 
quantum registers (cQPL and LanQ). 

• Optimisation of memory management (QCL). 

Both reasons are important from the simulations point 
of view, since they facilitate writing of correct and op- 
timised quantum programmes. However, these features 
do not provide a mechanism for developing new quantum 
algorithms or protocols. 
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